Exploiting Quench Dynamics in Spin Chains for Distant Entanglement and Quantum 

Communication 
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We suggest a method of entangling significantly the distant ends of a spin chain using minimal 
control. This entanglement between distant individual spins is brought about solely by exploiting 
the dynamics of an initial mixed state with Neel order if the lattice features nearest-neighbor XXZ 
interaction. There is no need to control single spins or to have engineered couplings or to pulse glob- 
ally. The method only requires an initial non-adiabatic switch (a quench) between two Hamiltonians 
followed by evolution under the second Hamiltonian. The scheme is robust to randomness of the 
couplings as well as the finiteness of an appropriate quench and could potentially be implemented 
in various experimental setups, ranging from atoms in optical lattices to Josephson junction arrays. 



The objective of quantum communication is to transfer 
a quantum state \<j>) from a sender (Alice) to a receiver 
(Bob) as accurately as possible. For this, Alice can sim- 
ply encode \<f>) on a carrier and send it. Alternatively, 
she and Bob can use teleportation for which they 
need to share a pair of particles in an entangled state 
= ^(U, T) + IT, I))- with appropriate opera- 

tions and classical communications, enables the noiseless 
transmission of a state \<j>) from Alice to Bob. Photons 
are ideal for long distance quantum communication. 

This paper is based on a manifestation of quantum 
communication which aims to connect distinct parts, i.e. 
registers, of a quantum computer. The idea is to use 
a many-body system (with permanent interactions) in- 
corporating the sender, the receiver, and the mediating 
channel all together. All components are then stationary, 
and could be made of the same units. This setup does not 
require the interfacing of stationary qubits with photons. 
This could considerably reduce the complexity of inter- 
connects between registers of a quantum computer which 
typically need to enable entanglement sharing over short 
distances i.e. several lattice sites. For the same reason, 
the physical movement of ions is considered seriously for 
communication between ion trap quantum registers 0]. 

Fig. [1] pictures our scenario, where Alice and Bob are 
situated at opposite ends of a one dimensional (ID) lat- 
tice of perpetually interacting spin-i particles. We sug- 
gest a scheme which allows the establishment of a strong 
entanglement between Alice's and Bob's spins (the re- 
motest spins of the lattice) without any requirement of 
local control for initializing the chain or for the sub- 
sequent dynamics. Neither the repeated switchings of 
any fields (local or global) nor any specially engineered 
couplings are required. This is close to the scenario of 
much work on state transfer through spin chains (e.g. 
Refs.0,0,[l,@,0|)» but here the lattice interactions "gen- 
erate" entanglement, as opposed to just "transferring" 
it. State transfer can itself be modified to yield entan- 
glement generation schemes (eg. see Q), but without a 
price (e.g. engineered couplings and local preparation), 
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FIG. 1: (Color online) Schematic of our proposal of entan- 
gling distant spins. Alice and Bob are at opposite ends of the 
chain. First, the spin chain is initialized in the antiferromag- 
netic Ising groundstate. By a non-adiabatic switching to an 
XY interaction and subsequent time evolution quantum cor- 
relations (entanglement) is being established between Alice's 
and Bob's spins. 



the amount of entanglement will be very small. Even 
without such price, our current mechanism provides a 
very high entanglement. 

In our scheme, first the lattice of strongly interacting 
spins is cooled to its ground state. Then, upon instantly 
changing a global parameter in its Hamiltonian (i.e., per- 
forming a quench 

[HEEQIH), the 

pair of edge 

spins evolve to a highly entangled mixed state. This state 
is special in the sense that entanglement purification [lt| , 
a process that Alice and Bob can use to convert, by lo- 
cal actions, a few (say n) copies of the state to m < n 
pure \ip + ) states, can be quite efficient. These could 
then be used to teleport any state from Alice to Bob. 
Though the scheme of this paper has a qualitative sim- 
ilarity with entangling the ends of a chain of uncoupled 
systems by a sudden switching of interactions [l(| , it 
yields a much higher entanglement scaling better with 
the length of the chain. Also, as we will show, we actu- 
ally start from a mixed initial state despite the cooling - 
so the production of a high entanglement between ends 
is quite interesting. 

Quenches have been actively studied in condensed mat- 
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ter, usually in the thermodynamic limit and were 
recently shown shown to make local subsystems relax to- 
wards maximum entropy In contrast, here we inves- 
tigate whether a quench, when performed on a finite sys- 
tem, can produce a "long-range" entanglement between 
individual distant spins. Recently, the development of a 
qualitatively different form of entanglement, namely that 
between complementary blocks of a s pin chain was shown 
to be possible after a quench 13, 14, [HI]. This phe- 
nomenon was attributed to a heuristic picture in which 
several sources in a lattice emit oppositely moving pairs 
of entangled quasiparticles which reach distinct blocks 
to entangle them. However, a high entanglement be- 
tween two individual spins is not guaranteed by the above 
picture - indeed for certain quenches in the quantum 
Ising model, where the above process entangles blocks 
[l3j |. there is hardly any entanglement between the end 
spins. A more intricate process involving reflections of 
wavepackets at boundaries, and a cumulative contribu- 
tion of the same entangled state from different sources 
to the end spins, can, for the quench we examine, create 
a very high entanglement between end spins at a special 
time. This is both directly relevant to quantum com- 
munications and measurable, as opposed to the entangle- 
ment of blocks. 

Using a global quench to create entanglement at a dis- 
tance is highly attractive because of the non-requirement 
of controlling single spins. Other low control ways of cre- 
ating sizeable entanglement between distant spins include 
exploiting certain spin chains whose end spins are cou- 
pled weakly to the remaining ones (l7l . [la ]. This is a 
rare engineered case, as generically the ground states of 
spin chain models are known to exhibit only very short 
ranged entanglement (HI- Another approach is based 
on localizable entanglement [5(| where one entangles the 
end spins of a bulk system by individual measurements 
on the other spins, which can be challenging. 

Returning now to the problem of entangling distant 
spins through a quench, consider a chain of N spin-^ 
systems with nearest-neighbor XXZ interaction 
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H =J2 yfe + ^ + A*!) (1) 
fe=i 

where the parameters Jk and A denote the coupling 
strengths at bond k and the anisotropy respectively, and 
CT fci a k an d a k denote the Pauli operators acting on the 
spin at site k. We assume Jk — J > (homogeneous 
anisotropic coupling) in most parts of this work, unless 
we study disorder in later parts of this paper: To that end 
the couping is modulated as Jk = J(l + Sk) with nor- 
mally distributed random numbers 5k with mean zero 
and standard deviation S. Eq. ([I]), corresponds to the 
Ising-Model for A — > oo, and the isotropic XY, or XX 
model for A = 0. As [H,S Z ] = 0, with S z = ^f =1 a%, 
the total z-magnetization is a constant of motion. 



We first formulate the analytic case of time-evolution 
of the Ising ground state, under action of the XX Hamil- 
tonian. This corresponds to an instantaneous, i.e. ideal- 
ized quench in the anisotropy parameter Ai — > A2 with 
Ai — * 00 and A2 = thereby crossing critical value 
A = 1, which separates the NeelTsing-phase from the 
XY-phase. Later on, in a purely numerical study we will 
allow for 1 < Ai < 00 and < A 2 < 1 . For A > J the 
Ising groundstate gets arbitrarily close to the ideal Neel 
state, which is twofold degenerate in the absence of an 
external field. These ideal Neel states arise from the per- 
fectly polarized state |JJ-at) upon flipping every other spin: 
|M> = Ui.ta.is,—} and \Af 2 ) = |Ti, I2, Ts, • • •> ■ Note 
that these two states turn into each other by a spin flip at 

each place, i.e. |A/i) = (IlfcLi a k)\^) an< ^ vice versa - I n 
an experiment, the initial preparation of the Neel-Ising- 
groundstate, will yield, at low enough temperatures, an 
equal mixture of both Neel orders, and negligible admix- 
ture of higher energy eigenstates. We adopt the notion 
of thermal ground state from [181 ] for 



Po = 



^(IMXMI 



W2W2 



(2) 



which exhibits the same symmetries as the Ising Hamilto- 
nian H (A — > 00), as opposed to each individual, degen- 
erate ground state of the antiferromagnetic Ising-chain. 

After diagonalization following a Jordan- Wigner trans- 
formation, we are provided the explicit time depen- 
dence of the local Fermi creation operators 



,t 



(n 



fe-i 
1=1 



a k , with a k 



l( a k + * °ic)' which reads 
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fk,i(t) 



N - 



(3) 

- sin(g m k) sin(q m I) e~ lEm * , (4) 



N 



7V+1 ' 



with E m — 2 J cos(g m ) and q„ 

We now investigate the dynamics of the long range 
nonclassical correlations, i.e. the entanglement between 
the first and the last spin of the chain, provided the sys- 
tem is initially prepared in the global state ((2|) and as- 
suming A = for the subsequent time evolution. The 
state (reduced density operator) pi : jv of spins at sites 1 
and N of the chain in the {|f |), \\ |), || |), || I)} basis 
has only the following non-zero elements 

(TT \pi,n\ TT) = (II \P!,n\ II) - a, (U \ P i.n\ TI> = 

(IT \pi,n\ IT) = 6, (TI \pi,n\ IT) = (IT \pi,n\ TI) = c, 

(5) 

which arc entirely expressible in terms of second moments 
of the Fermi operators ([3]). From a = (a^ cr^er^) and 
c = {(TiOjj), where (•••} = Tr{p(£) ••■} and pit) — 
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e lHt poe lHt , we find 



a = (cjci)i {J N c N )i - (cjcjv)i (cjvCi)i 
-2(( c ! c i)i + (c]vc w }i - 1) 

1 



(0) 



c = ^((-lr^c^ + ac.). 

Here, (--^i denotes the expectation value in the 
Schrodinger picture with respect to the initial state |A/i), 
which facilitates factoring higher moments according to 
Wick's theorem, and M is the conserved number of spin 
up states in the dynamical state \Afi(t)), i.e. M — TV/2 
for even N and M = (N - l)/2 for odd N . The re- 
maining element b is determined by the trace constraint 
2a + 2b = 1 . Analogous formulas for ground and ther- 
mal states of the XX spin chain were reported in [l9j |. 
The second moments occurring in ([6]) can be conveniently 
evaluated in the Heisenberg picture, where 

(#)<*(*)> = E (4(o)q(o)) . (7) 

As (4(0)q(0))i = S ki i8 k>2m , (m = 1, 2, ... , M), we get 



C =2<- 



M 

X E 

m= 1 



/jV,2m(*)/*,2m(*) + C ' C - 



(8) 



This expression underlines that quantum correlations be- 
tween the end spins are a superposition of contributions 
from spin flips originating from every other site 2m. 

We now examine whether a given supply of systems, all 
described by the same mixed state from Eq. © , can 
be converted to a smaller subset of states through 
local actions by Alice and Bob and classical messaging 
between them, for subsequent use in teleportation. A 
general criterion [l6[ for this procedure, called entangle- 
ment purification, to be possible for mixed state p of 
two qubits is expressed in terms of the fully entangled 
fraction f = max((e|p|e}). Here the maximum is taken 
with respect to all maximally entangled states {|e)}. The 
criterion reads / > i, and is adapted to our particular 
problem as max(6 + c,b — c) > i in view of Eq. §5§ under 
the trace constraint a + b = g. 

The function / will feature a first maximum some time 
Tmax after the instant of quench (Fig[2J left inset) and 
this is always found (by a numerical search for different 
TV) to scale as T max ~ (Fig. [2l right inset - in fact, 
Tmax is slightly lower than ^j). Interestingly, for chains 
with even TV and po as the initial state, pi,./v is always un- 
entangled (separable). This is because \JVi) contributes 
10 = a \ IT) + iP\ Tl) t° Pi.n (with real a and (3) while 
symmetry implies IA/2) to contribute an equal amount 
of |£) = a x ® a x \0 so that + is separable. 

Therefore, in Fig. [2] and henceforth, we consider odd TV. 

A high entanglement specifically between the end spins 
requires the reflection of wavepackets at the boundaries, 




FIG. 2: Fully entangled fraction / evaluated at time t = T max 
as a function of system size TV in a semilogarithmic scale. The 
data points for the analytic quench scenario (circles) exceed 
0.5 for chains up to TV = 241 spins. We find a good agree- 
ment of these data for TV > 25 to a function of the form 
g Rt {N) (x 1.42(9) N-" with v = 0.22(9) (dashed curve). For 
TV < 11 the figure is supplemented with data from numer- 
ical diagonalization for more general quench characteristics 
(Ai — > A2) : (00 — > 1) (triangles) and (3 — * 0) (squares) . Left 
inset: Typical time evolution of / in a quench (00 — > 0) for 
TV = 7 with first maximum at T max (vertical line). Disorder 
(S — 0,0.1,0.2,0.3 from top to bottom curves respectively) 
is seen not to affect the amount of entanglement seriously 
and leaves T max almost unchanged (average taken over 100 
realizations). Right inset : Linear scaling of T max with TV. 



the initial Neel order and an independent propagation of 
different fermions after quench. Consider |A/i) (with odd 
N) to have its central spin ||). Each |f) of |7Vi) (or each 
II) of IA/2)) is the location of a fermion which goes to a 
superposition of left and right moving wavepackets due to 
the dynamics. From Ref. 3, U| we know that such dy- 
namics alone suffices to entangle blocks, but entangling 



individual spins requires much more. T max is close to 
the time it takes a fermion to traverse half the chain. At 
Tmax a fermion initially, say, in the left half of the chain 
would have the peak of its left wavepacket reflected at the 
left end and traveled a distance d from that end, while 
the peak of its right wavepacket, yet to encounter an end, 
would be about the same distance d from the right end. 
This implies that for each fermion, the amplitudes of its 
presence at the two ends will be nearly equal at T max . 
A2 ~ (free model) can ensure an independent dynam- 
ics of each fermion. Both and ensure that each 
fermion crosses either an even or an odd number of oth- 
ers to reach either end, so that there is no net effect of 
exchange phases. Thus each fermion will independently 
contribute a state very close to to p\,N{T max ) (the 
nearer to an end a fermion is initially, the more its contri- 
bution deviates from \ip + )). The cumulative contribution 
of the same entangled state from all fermions at T max cre- 
ates a very high entanglement between 1 and N. If the 
different fermions contributed distinct entangled states 
(e.g. ||, t) + e^lt, I) with different 0's), the cumulative 
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entanglement could have been vanishingly small even if 
block entanglement would be high. 

We now provide an estimate for the supply of impure 
pairs required as input in order to produce one almost 
pure maximally entangled state, i.e. / out > 0.99. At 
time T max , our mixed state becomes p^jv — /|V' + }(V' + I + 
n ^ Q (lT,T)(T,TI + 11,1X1.11), and the so called recur- 
rence method of purification, described in detail in [l6j |. 
simplifies decisively (error correction, as usable for com- 
munication through disordered chains [2l| , could also be 
used). Starting from an ensemble of impure pairs with 
individual / = 0.544, which is the value for N = 151 
in Fig. [2] will require 5 iterations of the purification 
scheme on ~ 361 input pairs to achieve a single pair 
with / out = 0.996. In comparison, for a particular chain 
of N = 9 spins, which is also a representative number for 
possible experiments, we have initially / = 0.9117 and 
will need ~ 3 impure pairs to be purified into an almost 
perfect \ip + ) in a single iteration step until our thresh- 
old is exceeded by / out = 0.991. We have also studied 
numerically and plotted in Fig [2] more general quenches. 
The results suggest that rather than having a perfect 
Neel order in the initial state (Ai — > oo), it is more vital 
to propagate according to the free fermion Hamiltonian 
(A2 = 0) in order to achieve purifiable entanglement for 
larger N. Our scheme is also suprisingly robust towards 
the type of disorder under consideration, i.e. random 
couplings. As shown in the left inset of FigJ^l for ran- 
dom offsets as large as an average ten per cent of the 
coupling strength J the amplitude of the fully entangled 
fraction at T max as well as T max itself remain virtually 
unaffected. 

The verifiability of our results is within experimen- 
tal reach. Promising realizations of XXZ spin chains 
are Josephson junction arrays [221] (where one can tune 
J and A by varying voltages [23| or magnetic flux [24|). 
Spin-spin interactions in optical lattices [25j have recently 
been demonstrated [26| . These experiments also provide 
a global method of preparing the initial state: pairwise 
spin triplets \ip + ) are first prepared in the double wells 
of a superlattice, a magnetic field gradient is then ap- 
plied to attain a |f, J.) state in each double well, and 
finally, the long lattice instantaneously ramped down to 
create a Neel state. Ensembles of finite lattices should be 
simulable with superlattices and measurements of spin 
states are possible 26]. We have ignored decoherence 



for the moment in view of the relevant realizations (op- 
tical lattices are an ideal arena to study quenches in ab- 
sence of environments [27j |. and small N Josephson junc- 
tion arrays do not have significant decoherence over our 
(?max < 7-) time-scales [24|). Trapped ions are also 
known to simulate spin chain Hamiltonians and their 
non-adiabatic changes [28]. In future, it would be in- 
teresting to study the influence of a finite quench rate 
on the entanglement between the end spins. We started 
from a Hamiltonian whose eigenstates are product (un- 



entangled) states when written in terms of a local ba- 
sis i.e., the spin states at given sites of a lattice. We 
then quenched to a nearly free Hamiltonian which has 
momenta-like eigenstates. Thus a similar behavior may 
show in other models, e.g. Hubbard Hamiltonians when 
one quenches from an appropriate product of localized 
states to a free (such as the superfluid) parameter regime. 
Subsequent to the the work reported here, similar entan- 
glement dynamics has also been found by other authors 
[301 ] in driven ferromagnetic chains, which maps to our 
case in certain limits. 

Summary:- Non-equilibrium dynamics of a spin chain 
after a quench 11, 12, UJ, UA UM is an area of intense 
activity. The scaling of distant spin-spin correlations af- 
ter quenches is topical 13 ll |32j with interesting instances 
of polynomial scaling [3l|. Much harder is to make the 
truly "quantum" correlations (or entanglement) between 
individual distant spins "substantial" and scale polyno- 
mially, which we have shown to be possible through a 
certain global quench. This is important because the 
entanglement between individual spins in complex sys- 
tems, a quantity of wide interest is easily measurable 
through correlation functions, but is notoriously short 
ranged [18|. Our scaling of entanglement (~ N~°- 22 ) 
seems significantly superior to that in the statics or dy- 
namics of any spin model studied so far (except for spe- 
cialized couplings [13] )• Additionally, that we start from 
a mixed state and obtain so much entanglement between 
distant spins that it is efficiently purifiable to \ip + ), at 
least for small (N ~ 9) chains, and usable for linking 
solid-state based quantum registers is interesting. Tests 
of our results, e.g., in finite chains in optical superlattices 
[2^], seem feasible. 
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